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ABSTRACT 

With the XMM-Newton and Chandra observatories, we detected two ex¬ 
tremely embedded X-ray sources in the R Corona Australis (R CrA) star forming 
core, near IRS 7. These sources, designated as Xe and Xw, have X-ray absorption 
columns of ~3xl0^^ cm“^ equivalent to Ay ~180™. They are associated with the 
VLA centimeter radio sources lOE and lOW, respectively. Xw is the counterpart 
of the near-infrared source IRS 7, whereas Xe has no A-band counterpart above 
19.4™. This indicates that Xe is younger than typical Class I protostars, probably 
a Class 0 protostar or in an intermediate phase between Class 0 and Class I. The 
X-ray luminosity of Xe varied between 29< log Lx <31.2 ergs s“^ on timescales 
of 3-30 months. Xe also showed a monotonic increase in X-ray brightness by a 
factor of two in 30 ksec during an XMM-Newton observation. The XMM-Newton 
spectra indicate emission from a hot plasma with kT ~3-4 keV and also show 
fluorescent emission from cold iron. Though the X-ray spectrum from Xe is sim¬ 
ilar to flare spectra from Class I protostars in luminosity and temperature, the 
light curve does not resemble the lightcurves of magnetically generated X-ray 
flares because the variability timescale of Xe is too long and because variations 
in X-ray count rate were not accompanied by variations in spectral hardness. 
The short-term variation of Xe may be caused by the partial blocking of the 
X-ray plasma, while the month-long flux enhancement may be driven by mass 
accretion. 

Subject headings: stars: activity — stars: magnetic helds — stars: pre-main- 
sequence — stars: rotation — X-rays: stars 


1. Introduction 

Low-mass protostars are divided into two classes. Class 0 and Class I, according to 
their infrared (IR) and radio spectral energy distribution (SED). This classification gener¬ 
ally traces their evolutionary status. Class 0 objects are thought to be young {t ~10^ yr) 
protostars which mainly emit in the far-IR and submillimeter wavelengths with blackbody 
temperatures of <30 K (Andre et al. 1993). They are believed to be accreting mass dy¬ 
namically from their huge circumstellar envelopes. Class I objects are believed to be older 
protostars {t yr) at the end of the mass-accretion phase and emit in the near-IR at 

temperatures of T ~3000—5000 K. 



- 3 - 

Protostellar cores are generally hidden inside enormous gas envelopes. Hard X-rays can 
penetrate the thick molecular clouds, and hard X-ray observations have revealed high energy 
activity associated with Class I objects (Koyama et ah 1996; Grosso et ah 1997; Imanishi 
et ah 2001). The observed X-ray emission exhibits occasional rapid outbursts reminiscent 
of solar flares though current star formation theories do not predict solar-type magnetic 
dynamos in very young stars. Montmerle et ah (2000) proposed an alternative dynamo 
mechanism, in which fossil magnetic fields link the protostellar core with its circumstellar 
disk reconnect. 

Hard X-ray emission from other embedded sources was reported in the OMC-2/3 cloud 
(Tsuboi et al. 2001). The two detected sources show Class 0 characteristics: huge absorption 
{Nn ~ 1 — 3 X 10^^cm“^), no near-IR counterparts, and associations with millimeter radio 
clumps. However, follow-up radio and near-IR observations by Tsujimoto et al. (2004) did 
not unambiguously classify them as Class 0 protostars. In particular, one of these sources 
correlates with a centimeter radio source and a jet feature in the H 2 band, which indicates 
excitation by a jet from a nearby Class I protostar. Though Skinner et al. (2003) and Rho 
et al. (2004) have reported X-ray emission from other millimeter radio clumps in NGC 2024 
and the Trihd nebula, the photon statistics in those observations were too limited (<50 
photons per source) to identify their nature conclusively. To date no X-ray source has been 
clearly identihed with a bona-hde Class 0 object. 

The R Corona Australis (R CrA) cloud is a nearby star forming region [d ~170 pc, 
Knude & H0g 1998). Among many young stellar objects in the cloud, those in the IRS 7 
region have attracted particular interest as a site of ongoing star formation. The region 
contains double peaked, strong centimeter emission; an eastern peak is designated as lOE 
or IRS 7B, and a western peak is designated lOW or IRS 7A (Brown 1987; Feigelson et al. 
1998; Harju et al. 2001). It also contains two submillimeter peaks (van den Ancker 1999), 
multiple millimeter continuum peaks (Henning et al. 1994; Saraceno et al. 1996; Chini et al. 
2003; Choi & Tatematsu 2004), and a signature of strong bipolar outflows (Harju et al. 
1993; Anderson et al. 1997), but only one near-IR source, IRS 7 (Wilking et al. 1997, W97). 
These characteristics make the IRS 7 region a promising host of Class 0 sources. Koyama 
et al. (1996) detected hard X-ray emission and an intense flare from the IRS 7 region, which 
suggested the presence of a magnetically active protostar. 

In this paper, we report the X-ray detection of two extremely embedded sources in the 
IRS 7 star forming core during XMM-Newton observations in 2003. These results, combined 
with the analysis of two Chandra observations carried out in 2000 and 2003, and follow¬ 
up near-IR observations with the University of Hawaii 88-inch Telescope (UH88), help to 
determine the nature of the detected X-ray sources. 



2. X-ray Observations and Absolnte Position Correction 


We analyzed X-ray data obtained with the XMM-Newton observatory on 2003 March 
28 and 2003 March 29 (hereafter ObsxMMi and ObsxMM 2 , see Table 1 for details). The 
XMM-Newton satellite (Jansen et ah 2001) is composed of three nested Wolter I type X-ray 
telescopes (Aschenbach et al. 2000) with Enropean Photon Imaging Camera (EPIC) CCD 
detectors in their focal planes (Striider et ah 2001; Tnrner et ah 2001). The observations were 
obtained with the EPIC pn and EPIC MOSl and EPIC MOS2 detectors, in fnll frame mode 
with the medinm filter. We pointed at similar sky positions in ObsxMMi and ObsxMM 2 ; the 
IRS 7 star forming core was at 6' off-axis, where the 90% photon radins is ~1'. We analyzed 
the X-ray data nsing the software package SAS ver. 5.4.1 and HEAsoft ver. 5.2. We first 
processed the data nsing the SAS scripts “emchain” and “epchain”, and screened ont high 
backgronnd periods nsing standard criteria, exclnding events when the entire chip connt rate 
of “pattern = 0” events above 10 keV was >0.35 cnts s“^ for MOS and >1 cnts s“^ for pn, 
and exclnding events close to hot pixels or ontside the held of view. Finally, we selected 
events with pattern %4 for spectral analysis of the EPIC pn and with pattern %12 for the 
EPIC MOS and the other analyses of the EPIC pn. 

To help sonrce identihcation and investigate long-term hnx variations, we analyzed 
two archival Chandra observations. Chandra has a single high-performance X-ray telescope 
which provides snb-arcsecond imaging (Weisskopf et al. 2002). The Chandra observations 
were obtained nsing the Advanced CCD Imaging Spectrometer detector using the Imaging 
array (ACIS-I) on 2000 October 7 and 2003 June 26 (hereafter Obscxoi and Obscxo 2 , see 
Table 1). Garmire & Garmire (2002) presented the hrst results of Obscxoi- The IRS 7 
star forming region was at 2' and 0.' 2 off-axis in Obscxoi and Obscxo 2 , respectively, where 
the 90% photon radius is ~1". We analyzed the Chandra data using the software packages 
CIAO ver. 2.3 and HEAsoft ver. 5.2. For Obscxoi, we reprocessed the level 1 event data 
with “acis_process_events” to compensate for degradation of the CCD spectral resolution 
by charge transfer inefficiency (CTI, Townsley et al. 2000). For Obscxo 2 , we used pipeline- 
processed level 2 event data which was already corrected for CTI degradation. Finally, we 
selected events from both observations with the standard grades 0, 2, 3, 4, and 6 in the 
ASCA event grade system for further analysis. 

For the correction of the absolute coordinates, we measured positions of X-ray bright 
sources on the combined EPIC pn plus MOS images between 0.3—10 keV and positions of 
sources in the ACIS-I image between 0.5—9 keV using available source detection packages 
(SAS: edetecLchain, CIAO; wavdetect, see also Table 1). We correlated these positions with 
near-IR source positions in the 2MASS All-Sky Point Source Catalog (PSC) obtained from 
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the NASA/IPAC Infrared Science Archive (GATOR)^, whose astrometric accuracy is good 
to about 0." 1^. Each X-ray source has only one 2MASS counterpart within 4", so that 
the source identihcation is very reliable. After adjusting the X-ray positions to the 2MASS 
positions, the deviation of the Chandra and XMM-Newton coordinates from the 2MASS 
coordinates is < 0." 6. 


3. X-ray Image 

Figure 1 shows a “true-color” X-ray image of the R CrA star forming region, which 
combines EPIC pn and MOS(l-|-2) data taken in both ObsxMMi and ObsxMM 2 - Red, green, 
and blue colors in the image represent soft (0.2—1 keV), medium (1—3 keV), and hard (3—9 
keV) bands, respectively. The IRS 7 region is colored in dark blue, which indicates the 
presence of hard X-ray sources, possibly suffering strong low-energy absorption. The top 
left panel of Figure 2 shows a detailed view of the IRS 7 region of the hard band image. 
A strong X-ray source, labelled Xe, was detected with edetecCchain in both ObsxMMi and 
ObsxMM 2 (Table 2). Though edetecCchain failed to detect a small peak on the tail of the psf 
to the north-west of Xe, our subsequent analysis using the ObsxMMi data conhrmed it as a 
weak source at ~14cr signihcance which we labelled Xw- The positional accuracy of both 
sources is hard to estimate because the sources are unresolved. In this paper, we tentatively 
put the errors at ~2", the half pixel size of the EPIC pn camera. If we rely on the result 
of edetecCchain, the absolute position accuracy of Xe is 1." 7 for ObsxMMi and 1." 2 for 
ObsxMM 2 at 90% conhdence. 

Chandra detected two weak X-ray sources in the IRS 7 region at above 4cr signih¬ 
cance in both Obscxoi and Obscxo 2 (fop right panel of Figure 2 for Obscxoi; Table 2). 
Though each source has less than 20 photons, the absolute source positions are determined 
within 1" thanks to the excellent spatial resolution of Chandra. These detected sources have 
corresponding XMM-Newton sources though, if we take the error circles derived using the 
edetecCchain results, the source positions of Xe between the Chandra and XMM-Newton ob¬ 
servations have a signihcant discrepancy of ~2". Hereafter, we refer to the Chandra detected 
positions for both sources. 

The bottom panels of Figure 2 compare the X-ray source positions with sources detected 
at other wavelengths. Xe is associated with the VLA centimeter source lOE (Brown 1987) 


^ http://irsa.ipac.caltech.edu/applications / Gator / 

^Sect. 2.2 of Explanatory Supplement to the 2MASS All-Sky Data Release: 

http: //www. ipac.caltech.edu/2mass / releases / allsky / doc / 
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and the sub-millimeter peak vdA 5 (van den Ancker 1999). Choi & Tatematsu (2004) 
also suggested a marginal millimeter source at the position of Xe- Xw is associated with 
centimeter (lOW in Brown 1987), millimeter (source 4 in Choi & Tatematsu 2004), and 
sub-millimeter sources, and the near-IR source, IRS 7 (see also section 5). 


4. X-ray Light Curves and Spectra 
4.1. East Source (Xe) 

We extracted source photons from the XMM-Newton observations of Xe from a 27." 5 
radius circle centered on the source position excluding the region around Xw, with back¬ 
ground taken from appropriate source free regions. We obtained ~2,000 net counts from 
ObsxMMi and ObsxMM 2 , enabling us to perform detailed timing and spectral analyses. For 
the Chandra data, we extracted events from a 1." 8 radius circle centered on the source 
position. The background level was negligibly low. Only 10—20 counts were extracted from 
the Chandra observation, which precludes detailed spectral analysis. 

The top panel of Figure 3 shows the background subtracted EPIC pn-l-MOS light curve 
of Xe in the 2—10 keV band. The first half of the light curve, corresponding to ObsxMMi, 
is mostly flat with some indications of a slight increase at the end. A constant provides an 
acceptable fif at above 90% conhdence (Table 3). The source was about four times brighter 
than the average count rate of ObsxMMi at the beginning of the second half, corresponding 
to ObsxMM 2 , and the count rate gradually increased by a factor of two. This part of the 
light curve can be acceptably fit at greater than 90% conhdence level by a linear increase 
with a slope of ~9.3 xl0“^ cuts s“^ day“^(Table 3). During both observations, the hardness 
ratio dehned as count rates in the 5—10 keV band over those in the 2—5 keV band remained 
unchanged, except for a minimum at 7.86 days. This means that time variation in count 
rate was not accompanied by any signihcant change of the spectral shape. On the other 
hand, the light curves show marginal dips on timescales of <1 ksec near 8.02 and 8.17 days. 
These dips are seen in both the pn and MOS light curves of Xe and other nearby sources, 
such as R CrA, IRS 5, and CrA 1, have no dips at those times. These facts do not support 
an instrumental origin. We note that these dips might suggest a tentative period of ~13.9 
ksec, with the combination of minimal dips at 6.56 and 7.85 days (the dip at 7.85 day was 
only covered with the MOS data and does not appear in the Figure 3). 

The EPIC spectra of Xe in ObsxMMi and ObsxMM 2 (Figure 4) show several similarities: 
signihcant emission up to ~10 keV; strong absorption below 2—3 keV; a broad line feature 
between 6—7 keV; marginal lines between 5—6 keV in the EPIC pn spectra (which may 



7 


be of instrumental or cosmic origin^). To investigate the 6—7 keV line feature, we fit the 
EPIC pn and MOS (1+2) spectra simultaneously with an absorbed power-law model with 
a Gaussian component. An acceptable £t above 90% confidence has a photon index of 3.0 
(2.5-3.4), a Gaussian line centroid of 6.60 (6.53-6.67) keV, and a Gaussian width of 0.15 
(0.079-0.28) keV, where the numbers in parentheses denote the 90% confidence range. The 
derived Gaussian width, equivalent to An ~7,000 km s“^ if produced by Doppler broadening, 
is unreasonably large for a stellar plasma (see also discussion 6.1). We therefore interpret 
the broad feature as a blend of iron lines from a hot plasma at 6.7 keV and a fluorescent 
iron line at 6.4 keV though the profile needs to be confirmed with deeper observations. 
We then fit the XMM-Newton spectra with XSPEC by an absorbed 1-temperature (IT) 
optically thin, thermal plasma model (wabs: Morrison & McCammon 1983; MeKaL; Mewe 
et ah 1995) with a Gaussian component with line center fixed at 6.4 keV. For either of the 
two observations this model yields acceptable fits at the 90 % confidence level (Model A 
and B in Table 4). Ah differs significantly between ObsxMMi and ObsxMM 2 , possibly due 
to the simplistic spectral model we assumed. Indeed, the soft emission below 3 keV was 
unchanged between ObsxMMi and ObsxMM 2 , suggesting an additional component along with 
the hard emission. In the single temperature fit to ObsxMMi, the soft emission is included as 
a part of the absorption slope whereas the fit to ObsxMM 2 determines Ah from the 3—5 keV 
slope and does not reflect the soft emission. We therefore refit the spectra of ObsxMMi and 
ObsxMM 2 simultaneously with an absorbed 2T model — IT for the variable hard component 
and IT for the constant soft component — with a Gaussian line at 6.4 keV. In this model, 
we tied the Ah of the hard components in ObsxMMi and ObsxMM 2 and tied the elemental 
abundances of all components. We allowed Ah of the soft and hard components to be fit 
independently because a model fit with a common Ah gives large Ah ~2.4x 10^^ cm“^ and 
hence an unrealistically large intrinsic log Lx ~ 35 ergs s“^ for the soft component. The 
model, again, successfully reproduced the spectra above the 90% confidence level (Model G 
in Table 4). The derived physical parameters of the hard component are at the higher end 
among those of Glass I protostars (e.g., see Imanishi et al. 2001, for comparison): large Ah 
~ 2.8xl0^^ cm“^, equivalent to Ay ~180™ (using the Ah —Ay relation by Imanishi et al. 
(2001)); plasma temperature of 3-4 keV; log Lx ~30.8 ergs s“^ in ObsxMMi, which further 
increased to ~31.2 ergs s“^ in ObsxMM 2 ; and a fluorescent iron line equivalent width (EW) 
of ~810 (240—1400) eV in ObsxMMi and ~250 (100—400) eV in ObsxMM 2 - Meanwhile, 
the metal abundance is ~0.2 (0.1—0.3) solar, which is typical of low-mass young stars (e.g. 
Favata et al. 2003). 


^Sect. 3.3.7.2 in XMM-Newton Users Handbook, 
http://xmm. vilspa.esa.es/external/xmm_user _support/documentation/uhb/index..html 



Since we have very few counts from the Chandra observations, we simply calculated a 
softness ratio, dehned as S / (H + S), where S and H are photon counts in the 0.5-3 keV 
and 3-9 keV bands, respectively. The ratios are 0.15±0.08 (S = 3, H = 17) in Obscxoi and 
0.38±0.13 (S = 5, H = 8) in Obscxo 2 where the errors show 1-sigma. The same softness 
ratios evaluated for the XMM-Newton best-£t models after adjusting to the Chandra ACIS-I 
response are 0.11 for ObsxMMi and 0.048 for ObsxMM 2 , suggesting that the spectra have been 
softer during the Chandra observations. This is consistent with the picture that the hard 
component further declined during the Chandra observations while the soft component was 
unchanged. 


4.2. West Source (Xw) 

We extracted XMM-Newton source events of Xw from an ellipse with axes 15" by 10." 
5 elongated toward the NNE direction, excluding a region with strong contamination from 
Xe, and selected the background from a symmetrical region near Xe. We did not use the 
ObsxMM 2 data, because the Xw region suffered strong contamination from Xe. For the 
Chandra data, we extracted source events from a 1." 8 radius circle. The background level 
was negligibly low. 

The XMM-Newton spectra were reproduced by an absorbed IT model with N-^ ~3.4x 10^^ 
cm“^, kT >1.6 keV, and log Tx ~ 30.5 ergs s“^ (left panel of Figure 5, Model D in Table 4). 
We added both spectra from Obscxoi and Obscxo 2 because their count rates are about the 
same. The summed spectrum can be reproduced by the best-fit XMM-Newton model just 
by changing its normalization (right panel of Figure 5, Model F in Table 4). 


5. Follow-up Observations in the Near-IR 

Infrared source catalogues in the R CrA region currently available have several shortcom¬ 
ings: shallow limiting magnitude (e.g. Kumu ~ 15™ for 2MASS, 16.5™ for W97); insufficient 
spatial resolution (e.g. 1" pix“^ for 2MASS and 0." 75 pix“^ for W97); and insufficient 
positional accuracy (e.g., ±1" for W97). We therefore analyzed two deep X-band images 
of the field obtained in August 1998 and August 2003 using UH88 and the near-IR imager 
QUIRC. During the observations, we used the f/10 secondary mirror, yielding a pixel scale of 
0." 1886 pix“^. The spatial resolution of our images is 0." 5 and 0." 8 in FWHM for the 1998 
and 2003 data respectively. Both images were aligned with an accuracy of 0." 3 by referring 
to the coordinates of R CrA in the 2MASS catalogue and a high resolution image (FWHM 
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~0." 14) obtained with the SUBARU telescope. The analysis was made with IRAF^. 

Both images show R CrA as a hlamentary and mildly extended reflection nebulae (Fig¬ 
ure 6 for the 1998 data), but the location of Xe shows no significant emission, though it 
does have some marginal enhancement. The flux upper-limit was measured by subtracting 
the extended nebular emission, which we estimated with a third-order polynomial surface 
fit. The 5-cr upper-limit of Xe within the 20 x 20 pixel box centered at Xe was ~19.4™ 
for the 1998 image and ~19™ for the 2003 image. Meanwhile, using the SUBARU image 
(Nedachi et al. in preparation), we measured the absolute position of IRS 7 at ( 0 : 2000 , <^ 2000 ) 
= (19^1^55.^ 34, —36° 57' 21." 69) with an accuracy of 0." 3. This is much more accurate 
than the earlier observations. With this new position, IRS 7 falls within the positional error 
circle of Xw (see the bottom right panel of Figure 2). 


6. Discussion 

6.1. The Nature and X-ray Emission Mechanism of Xe 

Harju et al. (2001) suggested that the radio counterpart of Xe, source lOE, might be a 
radio galaxy or Galactic microquasar. However, AGNs with the observed X-ray flux above 
2—8 xl0“^^ ergs cm“^ s“^ (2—10 keV) are found <10 degree in the sky (Ueda et al. 1998), 
and the probability to detect such an AGN in the IRS 7 star forming core (~10arcsec^) is 
extremely small (<10“^). Furthermore, AGNs do not generally show thermal iron K emission 
line and have rather flat spectral slopes (F <2) (Ueda et al. 1998). Similarly, Galactic black 
hole candidates also show flat power-law X-ray spectra with F = 1.5 — 2.1 though they 
sometimes show thermal spectra with kT ~0.5-1.5 keV (McGlintock & Remillard 2003). 
Taking its association to the star forming core into account, Xe is most likely a very young 
stellar object. 

Gompared to Glass I objects in the R GrA star forming region such as IRS 1, 2, and 5 
which have Ah ~2x 10^^ cm“^ and K <11™' (from our additional analysis of the ObsxMMi 
and ObsxMM 2 data; see also W97), Xe shows much larger X-ray absorption and much smaller 
near-IR luminosity. Xe is also associated with the strong submillimeter condensation vdA 5 
(van den Ancker 1999). All these suggest that Xe is much younger than typical Glass I 
objects, and that it is a Glass 0 object or an object in an intermediate phase between Glass 0 


"‘IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As¬ 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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and Class I. We note that extremely embedded sonrces in the OMC-2/3 clond have similarly 
high absorption colnmns of A^h (1 — 3) X 10^^ cm ^ (Tsnboi et al. 2001). 

Between the Chandra and XMM-Newton observations, Xe exhibited strong long-term 
X-ray variation by a factor of 10—100 on a timescale of a month (Fignre 7). In none of 
the observations did we detect obvions flare activity thongh ObsxMM 2 showed a marked flux 
increase. Active stars such as RS CVn and young stars in open clusters do not generally 
vary in X-rays more than a factor of 2—3 outside flares (Stern 1998). Less active stars such 
as the Sun exhibit strong X-ray variations by up to a factor of 100, coincident with their 
activity cycles (e.g. Favata et al. 2004), but, unlike Xe, the X-ray luminosity of such stars 
is typically less than 10^® ergs s“^ and the observed activity time scale is several years. One 
possibility is that the strong variability of Xe could indicate abrupt activity produced by an 
enhanced mass accretion episode similar to that recently attributed to the outburst of the 
star in McNeil’s nebula (Kastner et al. 2004). Indeed, the outburst increased the X-ray flux 
by a factor of 50, and the post-outburst X-ray luminosity of 10^^ ergs s“^ is comparable to 
the luminosity of Xe during ObsxMM2- 

The plasma temperature and X-ray luminosity of Xe exceed the typical quiescent X-ray 
activity of Class I protostars and are comparable to temperatures and luminosities of X-ray 
flares from Class I protostars (Imanishi et al. 2001; Shibata & Yokoyama 2002)^. X-ray 
flares from Class I protostars may be produced by reconnection in a magnetosphere which 
is twisted due to the core-disk differential rotation (Tsnboi et al. 2000; Montmerle et al. 

2000) . Perhaps a similar mechanism explains the X-ray emission from Xe during the XMM- 
Newton observations, though magnetic reconnection would have to occur throughout the 
XMM-Newton observations since no rapid X-ray variation was seen from Xe. 

The flux increase of a factor of two in ~30 ksec in ObsxMM 2 is unlike the types of 
variations seen in magnetically driven X-ray flares which are characterized by rapid (~ 
10 ksec) flux increases (e.g. Tsnboi et al. 1998, 2000; Stelzer et al. 2000; Imanishi et al. 

2001) . Favata et al. (2003) found a similar rise in X-ray brightness in the classical T-Tauri 
star XZ Tan, with a factor of 4 increase during 50 ksec. In this case, the brightening was 
accompanied by an Ah decrease and therefore Favata et al. (2003) interpreted it as an eclipse 
of the emitting region by the accretion stream. Because Xe did not show any signihcant 
hardness ratio variation the absorber would have to be uniformly dense, optically thick gas. 
Such a variation could be produced by an eclipse of the X-ray emitting region by an absorber 


^Imanishi et al. (2001) used the distance to the p Oph cloud of 165 pc instead of 120 pc derived from 
more reliable HIPPARCOS data (Knude & H^g 1998) for a comparison to earlier X-ray results of the p Oph 
field. Their X-ray luminosity should be divided by a factor of two to compare to our result. 



or emergence of the X-ray emitting region from behind the rim of the protostellar core as a 
result of stellar rotation. 

If the fluorescent iron line in the spectra is real, this is unusual because fluorescent iron 
lines have been rarely observed from pre-main-sequence stars. Even a few examples during 
strong flares from Class I protostars have EW< 150 eV (Imanishi et ah 2001). The large 
equivalent width of the fluorescent line from Xe (~250—800 eV) again suggests that the 
source is extremely embedded. When we simulate fluorescent iron line EWs, assuming solar 
abundance for the surrounding cold gas (Inoue 1985), an optically thick absorber should 
block the direct X-ray emission by ~60% for ObsxMMi and ~3% for ObsxMM 2 - This result 
is consistent with obscuration of the X-ray emission though the blocking factor in ObsxMM 2 
should be >30% to explain the observed flux increase in ObsxMM 2 - Interestingly, the intrinsic 
X-ray luminosity in ObsxMMi should be log Lx ~31.2 ergs s“^, which is comparable to Lx 
in ObsxMM2- 

Our thermal model £t requires a metal abundance of ~0.2 (0.1-0.3) solar. Though 
the derived abundance is dependent on the thermal model, and identifying emission lines in 
low resolution spectroscopy may be difficult especially around ~1 keV (e.g. Kastner et al. 
2002), the metal abundance of Xe is determined mostly from the iron K line region, which 
seems to show a real underabundance. X-ray emission from T Tauri stars similarly show 
sub-solar abundances (e.g., Favata et al. 2003). This might suggest that Xe has an X-ray 
emission mechanism similar to T Tauri stars. Such abundance anomalies could be produced 
by the FIP (first ionization potential) effect (e.g. Giidel et al. 2001). Unfortunately, the 
spectrum of Xe does not have enough counts to detect, unambiguously, emission lines from 
other elements such as argon and calcium to test the FIP effect. Another possible solution 
is that the continuum emission includes a non-thermal component, as proposed for p Oph 
SI (Hamaguchi et al. 2003). 

From the above discussion, we hypothesize that the X-ray plasma responsible for the 
hard component was produced at a mass accretion spot on the protostellar core. The X- 
ray emission is blocked when the spot would have been behind the protostellar core during 
ObsxMMi, and just appeared from behind the rim in ObsxMM 2 as a consequence of proto¬ 
stellar rotation (Figure 8). To be consistent with the observed XMM-Newton light curves, 
the rotational period of the proto-stellar core would need to be >2.8 days. This rotation 
speed is much slower than the break-up rotation speculated for Class 0 protostars from 
rotational periods of Class I protostars (e.g. Montmerle et al. 2000). 

The soft component was apparently constant and had much smaller Ah compared to 
the hard component. This may suggest that the component has no physical connection to, 
and exists far from, the hot component. One possible origin is that the soft component 
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is associated with another hidden protostar, but, though the A^h of the soft component is 
typical of Class I protostars, the iC-band magnitude of >19™ is much larger than those of 
Class I protostars in the R CrA cloud {K <11™). Another possible origin is that the X-ray 
plasma is heated by a collision of a steady jet or outflow from Xe with circumstellar gas, 
a mechanism thought to be associated with X-ray emission from HH2, L1551 IRS 5, and 
CMC 2/3 (Pravdo et ah 2001; Favata et ah 2003; Tsujimoto et ah 2004). Indeed, Xe is 
associated with a centimeter radio source as those systems are, but the plasma temperature 
and X-ray luminosity of Xe are very large compared to those sources, except for the source 
in CMC 2/3. Such a high plasma temperature requires an energetic jet with Vjet ~1,500 
km s“^. While low-mass young stars generally have slow outflow velocities (a few hundred 
km s“^), Marti et ah (1995) measured a large proper motion in the young stellar jets HH80- 
81 implying velocities up to 1,400 km s“^. Xe could be another example of a source with 
high speed outflow. 


6.2. What is Xw? 

Xw is a counterpart to the near-IR source IRS 7 as well as the centimeter radio source 
low. The column density of Xw (~2—6x10^^ cm“^) is much larger than that typical of 
Class I protostars, while the X-band magnitude of ~12.2™ is comparable to the brightness 
of Class I protostars in the R CrA cloud. Therefore, Xw might be a Class I source seen at a 
large inclination angle. On the other hand, IRS 7 had been suspected to be the counterpart 
of an X-ray flare source seen during an ASCA observation (Koyama et al. 1996). However, 
the Ah measured during the ASCA flare (~4xl0^^ cm“^) was signihcantly smaller than 
that of Xw, even considering that the ASCA spectrum is contaminated by emission from 
surrounding Class I protostars with lower Ah. If the ASCA flare source is really Xw, the 
absorption column must be variable. 


7. Summary &; Conclusion 

We discovered two extremely embedded X-ray sources at the positions of the strong VLA 
centimeter radio sources in the IRS 7 star forming core. Thanks to its vicinity to the Sun 
{d ~170 pc), the large effective area of XMM-Newton and an opportunity to catch an active 
phase, we obtained around ~2,000 photons from Xe, which is about 40 times better than 
other extremely embedded X-ray sources observed so far. The upper-limit to the A-band 
luminosity was restricted to ~1/150 of the extremely embedded protostars in the OMC-2/3 
cloud {d ~450 pc and completeness limit ~16™, Tsujimoto et al. 2003). This is therefore 
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a rare, high signal-to-noise X-ray detection of what appears to be a Class 0 protostar or a 
protostar in an intermediate phase between Class 0 and Class I. The combination of XMM- 
Newton and Chandra observations enabled us to study the X-ray emission mechanism in 
detail. The source showed a significant long-term X-ray variation by a factor of 10-100. The 
X-ray properties - kT ~3—4 keV and log Lx ~31 ergs s“^- are comparable to magnetically 
active X-ray flaring sources though the light curve of Xe does not suggest flaring activity. 
These results may suggest that the X-ray activity was enhanced by mass accretion. The 
X-ray flux during ObsxMM 2 increased monotonically by a factor of two in 30 ksec, and the 
X-ray spectra showed a strong iron fluorescent line with EW of 250-800 eV. These phenomena 
may be related to the partial blocking of the X-ray plasma or perhaps they are indicative of 
rotation of a hot spot on the proto-stellar core. 

The evolutionary status of Xe should be constrained by constructing the IR and radio 
SED. This requires high spatial resolution and high sensitivity. Spitzer, ALMA, and 8 m 
class ground based telescopes are well-suited to examine this mysterious X-ray protostar. On 
the other hand, long exposure X-ray observations would confirm the fluorescent line profile 
and test the suggested periodicity in the X-ray brightness. These studies would help to 
understand the physical nature of protostars at the earliest phase. 

We are grateful to M. Choi, T. Hunter, T. Yokoyama, K. Tatematsu, M. Tsujimoto, 
H. Murakami, K. Koyama, S. Yashiro, Y. Tsuboi, and the anonymous referee for useful 
comments. This work is performed while K.H. held awards by National Research Council 
Research Associateship Award at NASA/GSFC, and is supported by XMM-Newton US 
grant. 
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Table 1. Observation Logs 


Observation 

Observatory 

Seq. ID 

Date 

Exposure 

(ksec) 

^axis 

Shift (Aa,A^) 

ObsxMMl 

XMM-Newton 

146390101 

2003 Mar. 28 

14.9/21.5 

6' 

(0." 4, -1." 0) 

ObsxMM2 

XMM-Newton 

146390201 

2003 Mar. 29 

18.1/24.0 

6' 

(0." 2, -1." 3) 

Obscxoi 

Chandra 

200017 

2000 Oct. 7 

19.7 

2' 

(0." 03, 0." 23) 

Obscx 02 

Chandra 

200194 

2003 Jun. 26 

37.6 

0.' 2 

(0." 42, -0." 03) 


Note. — Exposure: EPIC pn/MOS for XMM-Newton. Aaxis- Off-axis angle of the IRS7 region. 
Objects used for the position correction: IRS2, IRS5, HBC677, CrAl, ISO-CrAl36, HHlOl IRS 4, ISO- 
CrAl37, TY CrA and HD 176386. IRSl and ISO-CrAl34 were also used for the position correction of the 
Chandra data. 


Table 2. Detected Sources 


Source 

Designation 

Observation 


(q:2000i <52000)^ 

(h m s, d ^ ^^) 

Net cnts 
(cnts) 

i^-band 

(mag) 

Counterpart 

Xe 

XMMU J190156.3-365726 

ObsxMMl 

19 

1 56.29, 

-36 57 26.5 

467.7 

>19.4 

lOE"' . vdA 5(?)° 



ObsxMM2 

19 

1 56.32, 

-36 57 26.4 

1611.3 




CXOU J190156.4-365728 

Obscxoi 

19 

1 56.43, 

-36 57 27.7 

19.6 





Obscx02 

19 

1 56.39, 

-36 57 28.3 

12.7 



Xw 

XMMU J190155.3-365721 

Ot)SxMMl&;2 

19 

1 55.26, 

-36 57 20.8 

114.o'’ 

12.2° 

low'* , Source 4* , IRS7, vdA 3(?)° 


CXOU J190155.3-365722 

Obscxoi 

19 

1 55.34, 

-36 57 21.6 

12.7 





Obscx02 

19 

1 55.32, 

-36 57 21.8 

22.4 




Note. — Net cnts: EPIC pn + MOS for XMM-Newton. 

^'Positional uncertaintiess are for XMM-Newton and for Chandra. 

^Net counts in ObsxMMl- 

‘^Reference of the magnitude of R1 (IRS 7) to Wilking et al. (1997). 
‘^Reference to Brown (1987) and Feigelson et al. (1998). 

^Reference to van den Ancker (1999). 

^Reference to Choi Tatematsu (2004). 


Table 3. Fitting Result of the Xe Light Curve 


Observation 

Binning 

(sec) 

Model 

Constant 
(10“^ cnts s“^) 

Linear 

(10“^ cnts s“^ day“^) 

X'^/d.o.f. (d.o.f.) 

ObsxMMl 

2,000 

constant 

0.93 


1.63 (13) 

ObsxMM 2 

1,000 

constant 

5.3 


1.96 (26) 


1,000 

constant H-linear 

3.7^ 

9.3 

1.21 (25) 


Count rate at TJD = 12727.8564 day. 











Table 4. Fitting Results of the Spectra 


Source 

Observation 

Model 

Comp. 

Nil 

{10^^cm“2) 

kT 

(keV) 

Abundance 

(solar) 

log E.M. 
(cm-3) 

Flux6,4fceV 

(10“® cnts cm~^ s“^) 

X^/d.o.f. (d.o.f) 

log Lx 
(ergs s“^) 

Xe 

ObsxMMl 

A(1T) 


13.2 (9.6-18.3) 

5.1 (2.8-10.4) 

0.5 (0.2-1.4) 

53.2 (53.0-53.6) 

2.4 

1.22 (21) 

30.4 


ObsxMM2 

B(1T) 


24.7 (21.6-28.3) 

4.4 (3.4-6.1) 

0.2 (0.1-0.3) 

54.1 (54.0-54.1) 

3.4 

0.99 (66) 

31.2 


ObsxMMl &;2 

C(2T) 

Vari 

28.1 (22.5-33.8) 

2.7 (1.8-4.2) 

0.2 (0.1-0.3) 

53.8 (53.5-54.1) 

2.8 (0.8-4.8) 

0.97 (86) 

30.8 




Var2 

=Vari 

4.0 (2.9-5.3) 

=Vari 

54.1 (54.0-54.4) 

3.4 (1.4-5.5) 


31.2 




Const. 

4.2 (1.4-11.4) 

2.3 (0.4-) 

=Vari 

52.5 (51.8-53.5) 



29.6 

Xw 

CbsxMMl 

D(1T) 


33.7 (19.1-60.0) 

4.7 (1.6-) 

0.3 (fix) 

53.3 (52.8-54.3) 


0.65 (8) 

30.5 


CbscX01&;2 

E(1T) 


33.7 (fix) 

4.7 (fix) 

0.3 (fix) 

52.9 (52.8-53.0) 


0.60 (3) 

30.0 


Note. — • variable hard component in ObsxMMl (Ot)sxMM2)- Const.: constant soft component. Lx- absorption corrected X-ray luminosity in the 0.5 — 10 keV band. Distance 

assumes d ~170 pc. 
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Fig. 1.— XMM-Newton “true-color” image of the R CrA star forming region (ObsxMMi + 
ObsxMM 2 )- The image is color coded to represent hard band (3—9 keV) to blue, medium 
band (1—3 keV) to green, and soft band (0.2—1 keV) to red. The dotted rectangle shows 
the IRS 7 region (held of view of Figure 2 except the bottom right panel). Class I protostar: 
IRS 1, IRS 2, IRS 5, Herbig Ae/Be star: R CrA, TY CrA, HD 176386, Weak-lined T-Tauri 
star: CrA 1. 
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R.A. {J20i)0\ 


Fig. 2.— Magnified images of the IRS 7 region. The top left panel: the 3—9 keV EPIC pn 
+ MOS image of XMM-Newton combining ObsxMMi and ObsxMM 2 , the top right panel; the 
0.5—9 keV image of Chandra in Obscxoi, the bottom left panel: the XMM-Newton image 
overlaying flux contours in sub-millimeter at 450 ^m taken by van den Ancker (1999), the 
bottom right panel: source positions with their error circles including sources in near-IR and 
radio wavelengths. Diamonds show XMM-Newton sources, crosses Chandra sources, squares 
centimeter radio sources (Bnumber. Brown 1987, FCWnumber. Feigelson et ah 1998, Source 
positions are referred to Feigelson et al. 1998), x millimeter radio sources {CTnumber. Choi 
& Tatematsu 2004), a triangle near-IR source (IRS 7). The gray circle in the bottom left 
panel shows the beam size of the sub-millimeter telescope. Error circles in the bottom right 
panel are drown in dot-bar line for X-ray sources, dotted line for radio sources, and solid 
line for the near-IR source IRS 7. 
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as “count rates in the 5—10 keV band / count rates in the 2—5 keV band” {bottom panel). 
Points with error bars in the top and bottom panels are made from summed data sets of 
EPIC pn plus MOS (1+2). The horizontal axis is truncated Julian day (TJD) — 12720. The 
arrows show timings of the dip feature. Bins in the top panel have 2 ksec for ObsxMMi and 
1 ksec for ObsxMM 2 , and bins in the bottom panel have 4 ksec for ObsxMMi and 2 ksec for 
ObSxMM2- 
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Fig. 4.— EPIC pn {left) and MOS (1+2) {right) spectra of Xe in ObsxMMi {black) and 
ObsxMM 2 {red). The solid lines show the best-fit model of the simnltaneons fitting with 
EPIC pn and MOS. Dotted lines show the soft component, barred line the hard component, 
and dot-bar lines the Ganssian component for the line at 6.4 keV. Each bottom panel shows 
the residnals of the fit- 


































































































































23 



^ 1 2 5 10 ^ 1 2 5 

chcamel energy (keV) channel energy (keV) 


Fig. 5.— Xw spectra in ObsxMMi {left, solid - pn, dot - MOS) and Obscxoi plus Obscxo 2 
{right). Both panels show data points and their best-£t models (Model D and E in Table 4 
from left). 
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Fig. 6.— K-band image of the IRS 7 region obtained in 1998. The Xe and Xw positions 
measure with Chandra are at the left and right crosses, respectively. 
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Fig. 7.— Observed X-ray flux between 0.5-10 keV in long time scale {left Xe, right. Xw)- 
The vertical narrow bars on the data points show photon statistical error at 90% confldence 
level. The thick bar shows the variable range of Xe during ObsxMM 2 - 
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Fig. 8.— Possible geometry of the proto-stellar core, jet, and observer . The wavy lines show 
X-ray emission (solid black: direct X-rays, grey: fluorescent and scattering X-rays, dotted 
black: direct X-rays are partially covered). 



